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We examine the critical behavior near the metal-insulator transition (MIT) in the one-dimensional 
extended Hubbard model with the on-site and the nearest-neighbor interactions U and V at quarter 
filling using a combined method of the numerical diagonalization and the renormalization group 
(RG). The Luttinger- liquid parameter K p is calculated with the exact diagonalization for finite size 
systems and is substituted into the RG equation as an initial condition to obtain K p in the infinite 
size system. This approach also yields the charge gap A in the insulating state near the MIT. The 
results agree very well with the available exact results for U = oo even in the critical regime of the 
MIT where the characteristic energy becomes exponentially small and the usual finite size scaling is 
not applicable. When the system approaches the MIT critical point V — » V c for a fixed U, K p and 
A behave as |lnA|~ 2 = ca(V/V c — 1) and (K p — |) 2 = ex(l — V/V c ), where the critical value V c 
and the coefficients ca and ck are functions of U. These critical properties, which are known to be 
exact for U = oo, are observed also for finite U case. We also observe the same critical behavior in 
the limit of the MIT critical point U — > U c when U is varied for a fixed V . 

PACS: 71.10.Fd, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

A number of theoretical studies have been made on 
the one-dimensional (ID) extended Hubbard model with 
the on-site interaction U and the nearest-neighbor inter- 
action V as a simple model for quasi-lD materials 1-4 . 
It has been reported that this model shows a rich phase 
diagram including the metal- insulator transition (MIT), 
the phase separation, the spin-gapped phase and the su- 
perconducting (SC) phase 5-7 . In particular, the MIT 
at quarter-filling has attracted much interest as it takes 
place for finite values of U and V in contrast to the MIT 
at half-filling where the system is insulating except for 
U = 0. Therefore, the MIT at quarter-filling is impor- 
tant as a typical example of the quantum phase transi- 
tion caused by the electron correlation, and have been 
extensively studied by many authors 7-13 . However, the 
critical properties of the MIT have not received so much 
attention as it is more difficult to investigate the system 
in the limit of the MIT where the characteristic energy 
becomes exponentially small. In this paper, we wish to 
study the critical behavior near the MIT which is a typ- 
ical example of the quantum critical phenomena caused 
by the electron correlation. 

The extended Hubbard model is given by the following 
Hamiltonian 

H = -*I> t Wi" + h.c) 

i,<7 



where cj CT stands for the creation operator of an electron 

with spin a at site i and = c icr Ci a . t represents the 
transfer energy between the nearest-neighbor sites and 
is set to be unity (t=l) in the present study. It is well 
known that this Hamiltonian (1) can be mapped on an 
XXZ quantum spin Hamiltonian in the limit U — > oo. 
The term of the nearest-neighbor interaction V corre- 
sponds to the Z-component of the antiferromagnetic cou- 
pling and the transfer energy t does the A"-component. 
When the ^-component is larger than the A"-component, 
the system has a " Ising" -like symmetry and an excitation 
gap exists. For the Hubbard model, this corresponds 
to the case with V > 2t where the charge gap is ex- 
actly obtained 14 . On the other hand, in the case with 
"XY"-like symmetry (V < 2t), the system is metallic 
and the Luttinger-liquid parameter K p is exactly given 
by cos^) = -V/2 15 . 

In the finite U case, exact results have not been ob- 
tained except for V = 0. In this case, the weak coupling 
renormalization group method (known as g-ology) and 
the exact (numerical) diagonalization (ED) method have 
been applied. The <?-ology yields the phase diagram of 
the ID extended Hubbard model analytically, but quan- 
titative validity is guaranteed only in the weak coupling 
regime 1 ' 2,12,13 . On the other hand, the numerical ap- 
proach is a useful method to examine properties of the 
model in the strong coupling regime 4-11 . In particular, 
the numerical diagonalization of a finite-size system has 
supplied us with reliable and important information 5-10 . 

However, it is difficult for purely numerical approaches 
to investigate the critical behavior near the MIT where 
the characteristic energy scale of the system becomes ex- 
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ponentially small. To overcome this difficulty, we have 
recently proposed a combined approach of the ED and 
the RG methods. 9 ' 10 . This approach enables us to ob- 
tain accurate results of the Luttinger-liquid parameter 
K p and the charge gap A near the MIT beyond the usual 
finite size scaling for the ED method. The obtained re- 
sults of K p and A have been compared with the avail- 
able exact results for U = oo and found to be in good 
agreement 10 . The phase diagram of the MIT at quarter- 
filling together with the contour map of the charge gap 
A has been obtained on the U-V plane 10 . However, the 
critical behavior near the MIT was not discussed in the 
previous work. Here we extensively apply this approach 
to the critical regime of the MIT to elucidate the critical 
behavior of K p and A in the limit of the MIT. 

II. LUTTINGER LIQUID AND RG METHOD 

First, we briefly discuss a general argument for ID- 
electron systems based on the bosonization theory 1-3 . 
According to this theory, the effective Hamiltonian can 
be separated into the charge and spin parts. So, we turn 
our attention to only the charge part and do not consider 
the spin part in this work. In the low energy limit, the 
effective Hamiltonian of the charge part is given by 

where v p and K p are the charge velocity and the coupling 
parameter, respectively. The operator (f> p and the dual 
operator 6 P represent the phase fields of the charge part. 
g%± denotes the amplitude of the umklapp scattering and 
a is a short-distance cutoff. 

In the Luttinger liquid theory, some relations have 
been established as universal relations in one-dimensional 
models. 3 In the model which is isotropic in spin space, 
the critical exponents of various types of correlation func- 
tions are determined by a single parameter K p . It is pre- 
dicted that the SC correlation is dominant for K p > 1 

(the correlation function decays as ~ r Kp ), whereas 
the CDW or SDW correlations are dominant for K p < 1 
(the correlation functions decay as <~ r~( 1+Kp ^) in the 
Tomonaga-Luttinger liquid 3 . The critical exponent K p 
is related to the charge susceptibility % c and the Drudc 
weight D by 

Kp = \^ Xc Df , \ (3) 

with D = jj- 9 Q^i^ , where E (<j)) is the total energy of 
the ground state as a function of a magnetic flux N a cf) 
and N a is the system size 3 . Here, the magnetic flux is 
imposed by introducing the following gauge transforma- 
tion: cj nrJ — > e tm <t>/ N a c t^ f or an arbitrary site m. 




FIG. 1. The RG flow obtained from the numerical initial 
condition (solid lines) and that from the analytical one (bro- 
ken lines) for various U at V = 0. The filled squares are the 
numerical initial conditions with L\ — 8 and L2 = 12, and the 
open circles are the analytical ones. The filled diamonds on 
the K p axis are the exact results for U =3, 5, 8 and 00. Inset 
shows the RG flow of K p (€) (solid lines), the numerical result 
of K p (£) (filled circles) and the exact result (filled diamonds) 
for various V at U = 00. 

When the charge gap vanishes in the thermodynamic 
limit, the uniform charge susceptibility Xc is obtained 
from 



Ac E (N e + 2, N a ) + E (N e - 2, N a ) - 2E (N e ,N a ) ' 

(4) 

where E (N e , N a ) is the ground state energy of a system 
with N a sites and N e electrons. Here, the filling n is 
defined by n = N e /N a . 

We numerically diagonalize the Hamiltonian Eq. (1) 
up to 20 sites system using the standard Lanczos algo- 
rithm. Using the definitions of Eqs. (3) and (4), we 
calculate D and Xc from the ground state energy of the 
finite size system. To carry out a systematic calculation, 
we use the periodic boundary condition for N e = Am + 2 
and the antiperiodic boundary condition for N e = 4m, 
where A^ e is the total electron number and m is an integer. 
This choice of the boundary condition removes acciden- 
tal degeneracy so that the ground state might always be 
a singlet with zero momentum. 

At quarter filling, the 8kF umklapp scattering is crucial 
to understanding the MIT. The effect of the umklapp 
term is renormalized under the change of the cutoff a — > 
e l a. In this work, we adopt the Kehrein's formulation as 
the RG equations 16 ' 17 
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dK p {£) 
d£ 

dlogG(l) 
d£ 



G{tfK p {£f 

T(w p (i)-iy 

[2-8K p (£)}, 



(5) 
(6) 



where the scaling quantity £ is related to the cutoff a, 
T(x) is T-function and G(0) = g3±/2-Ka 2 v p . This formu- 
lation is an extension of the perturbative RG theory and 
allows us to estimate the charge gap together with K p 
in the infinite size system. To solve these equations con- 
cretely, we need an initial condition for the two values: 
K p (0) and G(0). Here, the value of the short-distance 
cutoff a is selected to a lattice constant of the system 
and set to be unity. Although the continuum field theory 
does not give this cutoff parameter, our choice is quite 
natural to apply this method to the lattice system. 

In the weak-coupling limit, analytic expressions for 
the initial condition have been obtained 13 . At quar- 
ter filling, v p , g 3 ± and K p (0) are given by {(2irvp + 
U + Wf - (U + W) 2 } 1 / 2 /2tt, (U - W)U 2 /(2irv F ) 2 
and {1 + (U + 4V)/(ttvf)}~ 1 / 2 , respectively, where vf 
is 2tsmkp- When we substitute these values into the 
RG equations as the initial condition, we find the in- 
sulating states for U > 8 at V = in contrast to the 
exact results which show the metallic states for all U at 
V = 0. This inconsistency suggests that the analytical 
initial condition is not applicable in the strong coupling 
regime. 

To find an adaptable initial condition to the strong 
coupling regime, we diagonalize a L-site system numeri- 
cally and calculate K p (£) by using the relation £ ~ InL. 
It is easy for the numerical calculation to obtain K p (£) 
as compared to G{£). Therefore, we calculate K p {£\) and 
K p {£.2) with L\- and L2-site systems instead of K p (£) and 
G{£) with a L-site system. To eliminate G{£) in the RG 
equations, we integrate the Eq. (6) and obtain, 



G(£) = G(£ a )J>- SK ^ M \ 



(7) 



where £o is a constant. Then the differential equation for 
K p {£) is written by 



d£ 



v(8K p (£) - 1) 



(8) 



When we set £q = £\ and use K p (£\) as the initial condi- 
tion for the above equation, we can obtain the solution 
numerically except the constant G{£\). By fitting the 
value of this solution at I = l 2 to _ftT p (f 2 ), we can deter- 
mine G{£\). Then, G(£) is immediately calculated from 
Eq. (7) and the solution of the RG equations is com- 
pletely obtained. 



III. DETAILED ANALYSIS OF K p 

In Fig. I, we show the RG flow obtained by solving 
the RG equations with the numerical and the analytical 




FIG. 2. Kp (oo) as a function of V at U = 5 10 and oo 
with the exact result at U — oo (the solid line). The broken 
lines are guides for the eye for the data of if£°(oo) . Inset 
shows Kp 2 (oo) as a function of for various V at U — oo 
(filled circles) together with the exact result (filled diamonds). 



initial conditions on the K p {£)-G{£) plane. Here, we set 
L\ = 8 and L 2 = 12 for the numerical initial condition. 
In the weak coupling regime with {/ < 5, the renormalized 
-Kp(oo) obtained from both the initial conditions agree 
with the Bethe ansatz result. On the other hand, in 
the strong coupling regime with U Z 5, there is a large 
discrepancy between K p (oo) with the analytical initial 
condition and the exact result. This is a striking contrast 
to the numerical initial condition which yields K p (oo) in 
excellent agreement with the exact result even in the limit 
U -> oo. 

In the inset of Fig. 1, we plot the RG flow of K p {£) 
as a function of L _1 for various V at U — oo together 
with K p (£) from the numerical diagonalization for sev- 
eral system sizes L and from the exact results for L — oo. 
The RG flow seems to connect smoothly the numerical 
results and the exact result. It indicates that the size 
dependence of K p {£) is well described by the RG equa- 
tions. For U = oo, our result is consistent with the pre- 
vious result from Emery and Noguera 18 . They solved 
the RG equations with the numerical initial conditions 
K p (£) and G{£) for a system size L, where G(£) is calcu- 
lated from the excited state energy. On the other hand, 
in our approach, we need only the value of K p {£) which 
is calculated from the ground state. Then, our approach 
can be easily extended to a complicated model such as 
the extended Hubbard model with finite U in contrast to 
the previous approach 18 which has been applied only for 
the infinite U case. 

In order to check the dependence of K p (oo) on the 
system sizes L\ and L 2 for the initial condition, we cal- 
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FIG. 3. The contour lines for K p x (oo) and for the charge 
gap on the U-V plane. The sold circles are K^(oo) from the 
RG equations with numerical initial conditions, and the open 
circles are the numerical diagonalization results (see in the 
text). The open squares represent the charge gap (Ref. 10), 
whose values are 0.25, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5 and 
5.5. The filled squares represent the phase boundary of the 
MIT where (oo) = 1/4. 

culate K p (oo) by using the three different sizes L 2 = 8, 
12 and 16 with Li = L 2 — 4. In Fig. 2, we show K p (oo) 
as a function of V for U = 5, 10 and 00 together with 
the exact result for U = 00. The value of K p {oo) is 
slightly dependent on L 2 . The inset in Fig. 2 shows 
the L 2 dependence of K p (oo) for various V at U — 00 
together with the corresponding exact result. Here, we 
assume that the size dependence of K p (oo) is given by 
K^ 2 (oo) ~ K™(oo) + Ci/Z| + c 2 /if, where a and c 2 
are constants, and K p x '(oo) is the L 2 — > 00 extrapolated 
value of K p {oo). We see that K p x '(oo) is very close to 
the exact result for U = 00. We may expect that the RG 
equations with numerical initial conditions give a reliable 
estimate for i^^°(oo) not only for the infinite U case but 
also for the finite U case where the exact result is not 
known so far. 

When the strength of V exceeds a critical value V c , 
K p (oo) is renormalized to the value of the strong cou- 
pling limit: K p — 1/4. This critical point corresponds 
to the MIT point of the system. In the infinite U case, 
we find V c ~ 1.93, 1.95 and 1.96 for L 2 = 12, 16 and 
20, respectively. Assuming the size dependence of V c to 
be cx I/L2, we obtain an L 2 — <■ 00 extrapolated value 
V c = 1.99. It agrees well with the exact value V c — 2 
for U = 00. The similar extrapolation yields the critical 
values of the MIT: V c ~ 3.45 for U = 5 and V c ~ 2.55 
for U = 10 as shown in Fig. 2 20 . The results are in good 



FIG. 4. I In A|" 2 as a function of (V/V c - 1) near the MIT 
at U — 00 for 1/2 = 12, 16 and 20 together with the exact 
result (solid line). Inset shows A as a function of V — V c - 

agreement with the phase boundary of the MIT in the 
previous works 5-10 . Then it confirms that the combined 
approach of the ED and the RG methods gives accurate 
results of K p even near the MIT. 

In Fig. 3, we show the phase diagram of the MIT on the 
U-V plane together with the contour lines for iGj°(oo) in 
the metallic region. We also plotted the contour lines 
for the charge gap in the insulating region which have 
already been reported in our previous paper 10 ' 19 . When 
V 3> U, the SC phase with K p > 1 appears. The char- 
acter of this phase has already been discussed in the pre- 
vious works 5-8 . Near the SC phase, K p (l 2 ) is larger 
than K p (£i) for available finite size systems and, then, 
we could not obtain the solution of the RG equations for 
these initial conditions. Because the umklapp scattering 
is canceled by the SC fluctuation, the RG equations may 
not be applicable in this region. Thus we estimate K p 
for V ^S> U directly by the ED method without the use 
of the RG method. 



IV. CRITICAL BEHAVIOR NEAR THE MIT 

Now we examine the critical behavior of the renormal- 
ized K p and the charge gap A near the MIT. In the per- 
turbative RG approach 8,21 , the asymptotic behavior of 
A is determined by the correlation length £ as A ~ v p /£, 
where In£ oc i -1 / 2 and t = (A— A c ) /A c . Here, A is a rele- 
vant parameter of a model indicating the deviation from 
the critical point A c . The critical behavior of K p is also 
obtained as (K p — 1/4) 2 cx t by analyzing the solution 
of the RG equations on the fixed line at G(oo) = with 
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FIG. 6. | In A|" 2 as a function of (V/V c - 1) near the MIT 
at U — 10 and (7 = 5 with L2 = 12 and 16 together with the 
straight lines whose slope are 0.20 for U = 10 and 0.22 for 
U = 5. Here, the factor of A is determined by fitting A to 
the numerical result at V — V c = 1. 



K p > 1/4. 21 

However, the derivation of £ in the perturbative RG 
approach is not clear because this approach leads the 
running coupling constants G(£) and K p (£) into diver- 
gence in the gapfull region with K p < 1/4. Further, the 
perturbative RG method fails to determine the explicit 
values of the prefactors of 't -1 / 2 ' and 'f in the formula- 
tion of ln£ and (K p — 1/4) 2 , respectively. 

On the other hand, in the combined approach of the 
ED and the Kehrein's RG methods, we can explicitly de- 
termine these factors without ambiguity. To avoid the 
difficulty of the divergence, Kchrein 16 ' 17 introduced a 
rcnormalizcd coupling constant G(£) constructed by the 
product of G(£) and the effective energy scale e~^ 2 ~ 8Kp \ 
In the limit £ — > 00, G{£) diverges in proportional to 
e £(2-8K p ) ^ gec ^)) ; w j 1 ji e rema i ns a finite value 
and is related to the charge gap as 

A = cv p G{oo), (9) 

where c is a factor of the order of unity. Here, the ex- 
plicit value c is estimated by fitting the RG result of 
G(oo) to the ED result of A and v p is calculated by the 
ED method. 22 The critical behavior of A including the 
explicit value of the prcfactor is determined by G(oo) and 
the critical behavior of K p is obtained straightforwardly 
from the running coupling constant K p (oo). 

A. V-dependence for U = 00 

First, we examine the critical behavior near the MIT 
for U = 00, where the MIT takes place when V is varied. 



In this case, we can test the reliability of our approach 
by comparing with the available exact results. Figure 
4 shows the critical behavior of the charge gap A cal- 
culated from the combined approach of the ED and the 
RG methods with L 2 = 12, 16 and 20 at U — 00, where 
I In A| — 2 is plotted as a function of (V/V c — 1) together 
with the exact result. In the inset in Fig. 4, A is plotted 
as a function of V — V c . Here, the factor c in Eq. (9) is 
determined by fitting the RG result of A to the numer- 
ical result from the ED method at V — V c = 1, where 
the system is away from the critical regime of the MIT 
and the ED method without the RG method can give an 
accurate result of A. 

As shown in Fig. 4, the critical behavior from our 
approach agrees very well with the exact result which is 
given by 

|lnV = CA (^" 1 )' (10) 

for < (V/V c - 1) < 1, where c A = ^- 14 . From the 
results shown in Fig. 4, we estimate the coefficient as 
c A - 0.087, 0.084 and 0.082 for L 2 =12, 16 and 20, re- 
spectively. Assuming the size dependence as oc 1/L|, we 
obtain an L2 — > 00 extrapolated value ca ~ 0.080 which 
is close to the exact result of ca = — 0.0821. Thus, in 
contrast to the perturbative RG method, our approach 
gives reliable and explicit estimates for the charge gap A 
with very small energy scale near the MIT. 

Fig. 5 shows (K p - 1/4) 2 as a function of (1 - V/V c ) 
at U = 00 calculated from our approach together with 
the exact result. In the critical regime, our result agrees 
very well with the exact result which is given by 
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FIG. 7. (JsTp - 1/4) 2 as a function of (1 - V/V c ) near the 
MIT at U = 10 and U = 5 for L 2 = 12 and 16 together with 
the straight lines whose slope are 0.019 at U = 10 and 0.032 
at U = 5. 

for < (1 - V/Vc) < 1, where c x = g^. From the 
results shown in Fig. 5, we estimate the coefficient as 
c K ~ 0.016, 0.015 and 0.014 for L 2 = 12, 16 and 20, 
respectively. Using 1/-L 2 , extrapolation, we obtain the 
coefficient in L 2 — > oo as ck ~ 0.013 which is again close 
to the exact result ck — g^y — 0.0127. This indicates 
that our approach gives an accurate estimate for K p as 
well as for A even near the MIT beyond the usual finite 
size scaling for the numerical diagonalization method. 

B. V-dependence for finite U 

Next, we examine the critical behavior of A and K p 
in the finite U case. In this case, there is no available 
exact result. In Fig. 6, we plot | In A|~ 2 as a function of 
(V/V c - 1) near the MIT for U = 10 and 5. We find that 
the critical behavior for both U = 10 and 5 is the same as 
that for U = oo given in Eq. (10) except the value of the 
coefficient ca- For U = 10, we estimate the coefficient 
as ca ~ 0.15 and 0.15 for L 2 = 12 and 16, respectively, 
which yield an L 2 — > oo extrapolated value ca ~ 0.15. 
For U = 5, the values of the coefficient are ca ~ 0.17 
and 0.18 for ^2=12 and 16, respectively, resulting in an 
L 2 — > oo extrapolated value ca ~ 0.19. 

Fig. 7 shows the critical behavior of K p for U = 10 and 
U = 5. We again find that the critical behavior for both 
U = 10 and 5 is the same as that for U = oo given in Eq. 
(11) except the value of the coefficient ck- We estimate 
the values of c K for U = 10 as c K - 0.022 and 0.022 



FIG. 8. | In A|~ 2 and (7f p -l/4) 2 as a function of \U/U C -1\ 
near the MIT at V = 5 for L2 = 12 and 16 together with the 
straight lines whose slope are 0.37 for | In A| -2 and 0.050 for 
(K p — 1/4) 2 . Here, the factor of A is determined by fitting A 
to the numerical result at U — U c = 1. 

for L2=12 and 16, respectively. For [7 = 5, the values 
of the coefficient are 0.028 and 0.029 for L 2 =12 and 16, 
respectively. These results yield L 2 — > 00 extrapolated 
values c K - 0.022 for U = 10 and c K - 0.030 for [7 = 5. 

When U decreases from U — 00, both ca and ck mono- 
tonically increase and become considerably large for a 
suitable value of U such as U — 10 and 5 as compared to 
the corresponding values of ca and ck for U = 00. 

C. [/-dependence for finite V 

In the large V regime (V > U), the MIT takes place 
at a critical value U c when U is varied for a fixed value 
of V as found in Fig. 3. Finally, we examine the critical 
behavior in this case. In Fig. 8, | lnA|~ 2 and (K p -1/A) 2 
are plotted as functions of \U/U C — 1| near the MIT for 
V = 5. We find that the critical properties as functions 
of U/U c are the same as those as functions of V/V c given 
in eqs. (10) and (11) except the values of the coefficients: 
I In A|- 2 - c' A (U/U c - 1) and (K p - ±) 2 = c' K (l - U/U c ) 
in the limit U — > U c , respectively. We estimate the co- 
efficients as c' A - 0.30 and c' K ~ 0.041 for L 2 = 12 
and c' A ~ 0.33 and c' K ~ 0.046 for L 2 = 16, which 
yield L 2 — > 00 extrapolated values as c' A ~ 0.37 and 
c' K ~ 0.050, respectively. Both of c' A and c' K are con- 
siderably larger than the corresponding values of ca and 

CK- 
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V. SUMMARY AND DISCUSSION 

In summary, we studied the critical behavior near the 
MIT in the one-dimensional extended Hubbard model at 
quarter filling by using the combined approach of the ED 
and the RG methods. In the large U regime (U > V), the 
MIT takes place at a critical value V c when V is varied 
for a fixed U, while, in the large V regime (V > U), it 
takes place at a critical value U c when U is varied for 
a fixed V. We examined the critical behavior near the 
MIT for both cases. 

In the large U regime, we observed the critical be- 
havior, | In A|- 2 = c A (V/V c - 1) and (K p - ±) 2 = 
ck(1 — V/V c ), where the critical value V c and the co- 
efficients ca and Ck are functions of U. For U = oo, the 
estimated values of V c , ca and ck agree well with the 
exact results. When U decreases from U = oo, all of V c , 
ca and ck monotonically increase. Both of ca and ck 
become considerably large for a suitable value of U such 
as U = 10 and 5 as compared to the corresponding values 
of ca and ck for U — oo. 

In the large V regime, we also observed the same criti- 
cal properties, | In A|^ 2 = c' A (U/U c - 1) and (K p -\f = 
c' K (l- U/U c ). Both of c' A and c' K are considerably larger 
than the corresponding values of ca and ck- For V ^> U, 
the SC phase with K p > 1 appears. Near the SC phase, 
it is difficult to obtain the solution for the RG equa- 
tion, because the umklapp scattering is canceled by the 
SC fluctuation. To examine the critical behavior near 
the MIT for V 3> U, we need an improved RG approach 
which includes both effects of the umklapp scattering and 
the SC fluctuation. 

We also obtained the phase diagram on the U—V plane 
and found that the phase boundary of the MIT and the 
contour lines of A and K p near the MIT smoothly con- 
nect between the large U regime and the large V regime. 
Although it is difficult to analyze the critical behavior for 
V ^> U, there is no qualitative difference in the critical 
behavior. These results suggest that the nature of the 
MIT is essentially unchanged on the U — V plane over 
the whole parameter regime including the large U and 
the large V regimes. 

In the limit V — oo, electrons are completely inhibited 
to occupy the nearest neighbor site of each other. In this 
case, some exact results have been obtained in the previ- 
ous works 5-8 : the ground state energy is always zero for 
U > U c (= 4) and the charge gap is given by A = \U—U C \. 
Then, the critical behavior of A for V — oo is completely 
different from that for finite V case obtained here. We 
think that there are two possibilities to explain this dif- 
ferences as follows: (1) The critical behavior near the 
MIT changes discontinuously at a finite V. (2) The criti- 
cal region smoothly shrinks with increasing V and finally 
becomes zero in the limit V = oo where the different 
critical behavior is observed. In any case, the critical 
behavior near the MIT for V U is interesting prob- 
lem as the competition between the SC fluctuation and 



the umklapp scattering becomes important, and will be 
studied in the future. 
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